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Abstract 

We use the general formalism of squeezed rotated states to calculate baryon asym- 
metry in the wake of inflation through parametric amplification. We base our analysis 
on a B and CP violating Lagrangian in an isotropically expanding universe. The B 
and CP violating terms originate from the coupling of complex fields with non-zero 
baryon number to a complex background inflaton field. We show that a differential 
amplification of particle and anti-particle modes gives rise to baryon asymmetry. 

1 Introduction 

Baryon asymmetry is one of the challenging problems in the standard model of cosmology 
PU |2] . It is characterized by the ratio rj = —, where n& is the number of baryons and n 7 
is the number of photons in the universe. The present value of the asymmetry in the uni- 
verse is f] ~ 10~ 10 . Three conditions postulated by Sakharov |3] are sufficient to guarantee 
baryon asymmetry. These are (i) Baryon number(B) violation (ii) charge and parity (CP) 
violation (iii) and the presence of nonequilibrium processes. These conditions are satisfied 
at very high temperatures, i.e. at Grand Unified Theory (GUT) scales. Thus the first the- 
ory of baryogenesis proposed was that of GUT baryogenesis [U E] . In GUT baryogenesis, 
baryon asymmetry is generated through the decay of baryon number violating bosons known 
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generically as X bosons. The general feature of GUT baryogenesis is that the characteristic 
energy scale for the processes with B-nonconservation is extremely high, close to the Planck 
scale, Mqut = 10 16 CreV, thus, it would be kinematically difficult to produce these Bosons 
in a thermal environment. At temperatures higher than the GUT scale, the rate of their 
production would be smaller than of the expansion of the universe and any generated baryon 
asymmetry is presumed to be wiped out by this expansion. If the temperature of the uni- 
verse could always be smaller than the GUT scale, the corresponding, thermally produced 
X-bosons would never have been abundant and their role in baryogenesis would be negligi- 
ble. Furthermore, together with B violating interactions that produce asymmetry in GUTS, 
there should be proton decay and neutron antineutron oscillations, which have, unfortu- 
nately, not been detected experimentally. Since then, many other scenarios of baryogenesis 
have been developed with baryon non- conservation at much lower energies. One such theory 
is electroweak baryogenesis where essential physical processes take place at an energy scale 
of around 100 GeV. The problem with electroweak baryogenesis is that the CP violation 
that generates the asymmetry is at least 12 orders in magnitude smaller than the observed 
asymmetry. One proposed solution to overcome these difficulties with baryogenesis models is 
to consider baryon asymmetry in the context of inflation. It is believed that an inflationary 
phase precedes the hot phase of cosmological evolution. In the inflationary phase, the energy 
density is dominated by the contribution from scalar field(s) called inflaton(s) A phase 
transition from an inflationary phase to a thermal phase may provide, in principle, the large 
temperatures needed for baryogenesis. During inflation, all matter and radiation are inflated 
away. Thus when inflation is over, there must be a reheating process, when the energy stored 
in the inflaton(s) is converted into thermal radiation. There are several possibilities for this 
reheating process. The inflaton may decay perturbatively into light particles, which thermal- 
ize eventually. There are also non-perturbative models of reheating known as "preheating " 
models. These are of two types. In one type of model the process of parametric resonance 
[HI E] may greatly enhance the production of particles ^U] • I n another, a rapid quench in the 
motion of the inflaton may lead to a spinodal instability and results in fast tachyonic preheat- 
ing [TI]. The preheating process is a far from equilibrium process forcing the universe into a 
non-equilibrium state after the end of inflation and before thermalization takes place. Thus 
it fulfills the third criterion for baryogenesis. It has been shown that some baryon-number 
violating processes similar to sphaleron transitions take place at preheating [T2*l IT3] , as well 
as tachyonic preheating ^3]. An interesting feature of the preheating scenario is that the 
production of particles with masses greater than the inflaton mass is possible. In particular, 
the non-thermally produced particles could include the baryon number violating gauge and 
Higgs bosons of grand unified theories (GUTs) or right-handed neutrinos, resurrecting the 
possibility of GUT-scale baryogenesis or leptogenesis . 

Preceding all of this, Papastamatiou and Parker proposed a simple model in which 
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the matter-antimatter asymmetry created by a baryon and CP violating Lagrangian in an 
isotropically expanding universe is found by a perturbative expansion. Taking an initial state 
in which no particles participating in the asymmetric interaction are present, they found that 
matter-antimatter asymmetry is produced during a stage when the radius of the universe is 
small with respect to its present value. In this work, we have adapted their work to consider 
the asymmetry that might be created in a more realistic universe that is initially inflating 
and then enters a rapid preheating phase either by parametric resonance or tachyonic pre- 
heating. Furthermore we use a more realistic source of CP violation, namely, a complex 
time varying background field which we can identify with a classical inflaton field , like in 
ref. [17J. Preheating is nothing but particle creation in the oscillating background, or one 
with spinodal instability (tachyonic preheating). Although particle creation or annihilation 
generally occurs in a time-dependent background, the number of created particles grows 
exponentially either when the background is periodic in time or when the coupling constant 
changes suddenly. This results in far from equilibrium particle production and there is an 
amplification of long wavelength modes. Therefore, if we add the ingredients of CP and B 
violation, there are chances for generating the matter-antimatter asymmetry of the universe 
in these eras. Recently, Ref. [T7| [TH] appeared in which the authors discuss the generation 
of the baryon asymmetry during reheating. In one of these ^B] baryon asymmetry was gen- 
erated by considering particle creation in the interaction of scalar fields with an oscillating 
inflaton background. In both these papers, baryon asymmetry was evaluated perturbatively, 
and the perturbation theory limited the produced asymmetry. In this paper, the first of a 
series, we extend the work in [T7| [TH] to include both expansion and non-perturbative effects 
and look for baryon asymmetry in the wake of inflation through parametric amplification. 
We calculate the effect of CP violation without resorting to perturbation theory. This allows 
us to compute non-perturbative out of equilibrium generation of baryon asymmetry. We do 
so by considering the Lagrangian of Papastamatiou and Parker in an expanding inflationary 
universe and treat particle production in the general formalism of squeezed rotated states. 

2 The Formalism 

The most general Lagrangian density with a baryon non- conserving interaction in an ex- 
panding inflationary cosmology has been given in Ej as 

+g»"dprfd u 7 } - {m% + C V R)V*V ~ V(v) ~ Kv 2 ^ + V*M)], (1) 

where, m^ f are masses of fields (p and ip. is a complex scalar field that carries a baryon 
number 1 and ip is a complex scalar field that carries a baryon number (which can be 
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treated as an antilepton field). The field 77 is a minimally coupled complex inflaton field, 
R is the Ricci scalar and £ is curvature coupling (£ = | gives the traditional conformal 
coupling of scalar fields to gravity). The sole purpose of the inflaton field 77 is to provide for 
inflation and therefore, it is treated as a background classical field while quantum excitations 
are considered for the other scalar fields. Our prime interest in this paper is to explore the 
effects of the inflaton field and the expansion of the Universe on the baryon assymetry. 
To this end, we consider the background gravitational field to be that of a standard FRW 
cosmology with an expansion parameter a(t). This Lagrangian can be thought of as an 
effective Lagrangian which contains terms that violate baryon number by means of a complex 
field which carries a baryon number interacting with the background inflaton field. This 
results in a net asymmetry between the produced particles and antiparticles represented by 
these complex fields, which gives rise to a net baryon asymmetry [T5] . It is this particle 
-antiparticle asymmetry that we now proceed to calculate. 

In order to see the salient features of baryogenesis, we simplify the Lagrangian density 
through the following assumptions. The inflaton field rj is a classical homogeneous time 
dependent complex background field, whose temporal evolution is determined by a Euler 
Lagrange equation of motion, with the inflaton potential V(rj) in the Friedman-Robertson- 
Walker metric ds 2 = dt 2 — a 2 (t)dx 2 , given by 

a dV 

77 + 3-77 + — = 0. 2 
a orj 

For this paper, we consider the minimal coupling to the background gravitational field = 
£v = £»? = 0)- The action is given by 

S = J d 3 xdta 3 (t)(±(<P*)(<i>) - (V<f ) (V0) + - ^(VV*)(VV0 

- m$*4>-m^*i>- AfaVV' + ^VW (3) 
Applying a Legendre transformation we can write the Hamiltonian as: 

H = J d^dta 3 [^tt;^ + ^ V0* V0 + ^vr;^ + ^ V^* W> + m|0*0 + 

+ A(t 7 2 ^V + 77*VV)] ) (4) 

where, is^, and ir^ are the canonical momenta of the <fi and ip fields. We canonically quantize 
this Hamiltonian. Since we are treating the inflaton field as a classical complex background 
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field we quantize only the 4> and ip field by introducing the mode expansions 

ij(x^) = J d~k[ate- ik - x + bfe lk -% 



= I dk[afe lk - x + bte~ ik - x } 
ke -ik-x + b W e ik. x y 



0(^) = J dk[a 
<i>*{xP) = J dk[afe lk - x + bt 



e ik ' x } 



here 



d 3 kdt 



dk = 



[(2tt) 3 2^; 



1 ! 
2 



and it is clear from the line element given above, the signature of the metric is (+,-,-, 
The quantized Hamiltonian is given as 

zj f j3; r^V \il> il> , u 4> \4> 4> i ^V 2 ® 3 ftp 6 . A?7* 2 a 3 u ^. 

y L a 3 fe k a 3 k k k k k kl 



where 



4 = i 2 + < 

fi - p + m 2 

„6 — +m V" 



here = - is the physical wave (co-wave) number of the mode, k is the comoving 
number. 

We define the following 




and write the Hamiltonian as 



k k a 3 



+& b * b t* + ^ b ffi + ( !!± 



+ 



a 



u k w k 



- 1 



w.1 



+ (16) 



To see the symmetries of this Hamiltonian, which will aid us in its diagonalization, we 
define the operators: 



JVi = a*a£ N 2 = a?al N 3 = b^'l N A = b?bl, 



(17) 



J + = a*a* J_ = a?a* J = -(N 1 -N 2 ), 
M + = bfbt M_ = b[X, Mo = \(N 3 - JV 4 ), 

satisfying the su(2) algebra 

[J + , J_] = — 2J , [J + , J ] = — J+, [J_, Jq] — 



(18) 



(19) 



and 



[M+, M_] = -2M , [M+, M ] = -M+, [M_, M ] = M_ 



= a*?b% K_ = b*_ k at K = l(N, + N 3 + 1), 



L + = 4*b%, L_ = b\al L = -(N 2 + N A + 1), 

satisfying the su(l,l) algebra 

[K+,K.] = 2K , [K + ,K ] = K+, [K.,K ] = -K_, 



(20) 



(21) 



(22) 
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[L + , L_] = 2L , [L + , L ] = L + , [L_, L ] = -L_. 
In terms of these generators we can write the Hamiltonian as: 



H = 



r 2uo. 



2uj 



a 3 



a;' 



col. 



(23) 



(24) 



The su(l,l) and su(2) symmetries of H are now manifest. We can digonalize this Hamiltonian 
by using the following unitary transformation 

H' = U\R 2 )U\R 1 )HU{R 1 )U{R 2 ), (25) 

where 

U(R 1 )U(R 2 ) = exp[9(J + e 2i ^ + J_e 2i ^)]exp[9{M + e 2 ^ + M_e 2i % (26) 
and operator U(R\) provides the well known transformation relations: 



U\R,) ( a a % ) U(R{) = 
while U(R 2 ) provides the relation 



cos(#) e 2i? sm(0) 
-e- 2i «sm(#) cos(#) 



cos(#) e 2t ^sin{9) 
-e~ 2l ^sin(9) cos(9) 



Bk 



Ck 



(27) 



(28) 



and the angle 9 is determined from the relation tan{9) = ^uj^/uj^ so that cos(9) = 
yjuj^jiuj^ + oity) and sin{9) = ^uj^/iu^ + u^). The inflaton field is r/ = \r]\e^ and 77* = 



1 77 1 e" 



After the rotations 1 the Hamiltonian assumes the form: 

3 (uj, + u;^) 



H' = U\R 2 )U\R 1 )HU{R 1 )U{R 2 



I 



d 6 k- 



[P(t)[AlA k + C^_ k ] 
+a(t)lBl k B_ k + D k Dl] 



a 3 



a 3 



w2 



1 



[Si^j + C?-*^]]. (29) 



1 This gives us the quantum optical analogy which we are exploiting in this work, as it is analogous to the 
effect of a lossless beam splitter, r = cos(9) and t = sin(6) are the reflection and transmission coefficients 
and the phase shift £ is between the reflected and transmitted fields. 
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This Hamiltonian has a su(l,l)x su(l,l) dynamical symmetry seen by defining: 

D 1+ = A\d k , D,_ = C^ k A k , D w = \{A{A k + d k C- k + 1), 

D 2+ = Bl k Dl D 2 _ = D k B_ k D 20 = \{Bl k B_ k + D\D k + 1), 

satisfying the su(l,l) algebras 

[D 1+ ,D 1 _] =2D 10 , [D 1+ ,D 10 ]=D 1+ , [D 1 _,D 10 ] = -D 1 _, 

[D 2+ ,D 2 _] = 2D 20 , [D 2+ ,D 20 ] = D 2+ , [D 2 _,D 20 ] = -D 2 _. 
Once again this symmetry is manifest through rewriting the Hamiltonian as 



(30) 
(31) 

(32) 
(33) 



/ 



CUT, 



■[\P(t)D 10 + a(t)D w ] 



[Di + + D 2 _] + ^ 



[D 2+ + D 1 ^)], 



where: 



2^6 



2a = 1- 
2/5 = 1 + 



A. j 77 J a 
2ui,puj^ 
A|^| 2 a 6 



(34) 

(35) 
(36) 



Since an su(l,l) symmetry implies diagonalizability through a squeezing (Bogolubov) trans- 
formation from the quantum optical analogies that we are applying, we see here that a 
product of two squeezing transformations will provide us with a diagonal Hamiltonian. The 
squeezing transformations we use are given by the following product of squeezing operators 

5(Ci)5(C 2 ) = exp[dD 1+ - £D 1 _]exp[{ 2 D 2+ - QD 2 _], (37) 

where, Ci = riexp[i^i] and ( 2 = r 2 exp[i^ 2 ] are the squeezing parameters. 
The operators 5 , (Ci)and S(( 2 ) provide the relevant Bogolubov transformations for the Hamil- 
tonian given in equation (29) in terms of the creation and annihilation operators A k , B k , C k 
and D k etc. 

A 9 (k,t)= f i 1 A k + u 1 Cl k , (38) 
Al(k,t)=riAl + v* 1 C- k , (39) 
B s (k,t)=fi 2 B- k + v 2 Dl (40) 
Bl(k,t)=^ 2 Bl k + u*D k , (41) 
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where ui = cosh(ri) = — M*L__ v _ e 171 s i n h(n) = e 111 2a ^' \ , U2 = cosh(r 2 ) = 

(4a(t)-l)~z (4a(t)-l)2 

— i and = e - * 72 sinh(r 2 ) = e~ n2 2/3 ^)~ 1 1 . 

(49(t)-l)3 (4/3(t)-l)2 

It is to be noted that obvious mode(k) dependence has been suppressed and will be made 
explicit when required. Thus the final diagonalized Hamiltonian is 

H f = J <Pk K 2 ^ } (4a(t)-l)* *M*ft «)+!]+ ( "^, faV) (49(*)-l)* 

(42) 

Let |0(t),0(t) > be the vacuum state of Hf(t). Then, the su(l,l) symmetry implies that its 
evolution is generated by the product of squeezing operators S((i) Sfa) given above. We 
can write this in terms of the generators of the su(l,l) x su(l,l) symmetry as: 

|0(t),0(t) >= e 1 ^¥^ Dl+ - Dl - )Ca{Da+ - D2 - ) \0(k),0(k) > (43) 

where \0(k), 0(k) >, vacuum state for H is defined as af. |0(fc),0(fc) >= 0, at\0(k),0(k) >= 0, 
bt k \0{k), O(Jfe) >= 0, and &^|0(ife), 0{k) >= 0. 

For a given Lie algebra, a Baker- Campbell-Hausdorff formula can be used to simplify the 
exponentials of operators in (jjSj) and is called disentangling in quantum optics j2Hj. The 
disentangling formula for su(l,l) is 

e Ci(Oi + -Oi-) e Ci(Oi + -Oi-)| (^) ) 0(A;) >= e 7lfcDl+ e r? i Dlo e 7 ^ Dl -e 72fcD2+ e r '2 D2( 'e 7 ^ D2 -|0(A;),0(A;) >(44) 

where, j lk = tanh(( lk ) and = 2ln(cosh((i k )) = -ln{\ - |7i fc | 2 ), i lk = 7* fc , 72/c = tanh(( 2 k) 
and 77' = 2ln(cosh((2k)) — —ln(l — |72fc| 2 ), l'i k = l 2k an d the modes are independent. This 
allows us to express the time evolution of the state as follows 

|0(t),0(t)> 



COsh(Cl(t))cOsh(C2(t)) 

]Jexp(tanh(( 1 (t))D l+ )exp(tanh(( 2 (t))D 2 +)\0(k),a(k) > . (45) 

k 

This is just a generalized su(l,l) squeezed state[26j. 

QFT in curved space times has taught us the importance of Bogolubov transformations and 
consequent population of one vacuum with particles of another vacuum at an earlier time. 
Here, we see this phenomenon with respect to both the squeezing (su(l,l)) symmetries. We 
see that the vacuum at times(t) is populated with particles and anti-particles with respect 
to vacuum at t=0. 

The number of particles and anti-particles can be calculated by the relationship between the 
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creation and annihilation operators of the initial quanta a%, b k , a\ and bf, before reheating 
with the final creation and annihilation operators A s and B s given by 

A s (k, t) = (pL lC os(9))at + (v lS in9e 2i Z)b^ k + (^sirafle 2 *)^ + (v lC os(9))b%, (46) 
B s (k, t) = {fi 2 cos{6))at k + {v 2 sin9e- 2 ^)b^ + {^sin9e-^)at k + (y 2 cos{0))b™ . (47) 

Note that it is the <fi field which carries the baryon number. The number of baryons at 
time t is then given by 



N B (t) = J2( A t(k,t)A s (k,t)) = Kil 2 , 

k k 

while the number of anti-baryons is given by 

k k 
Therefore we have a baryon asymmetry which is given by 

N B (t) - N^(t) = J2(Wki\ 2 ~ Wk 2 \ 2 ), 



(48) 



(49) 



(50) 
(51) 



where the k dependence is through u^, oj^ 

For A < 1, i.e the perturbative limit we recover the results of Papastamatiou and Parker, 
the asymmetry being zero to orderA, and non zero at order A 2 and above. 



Aa 6 M 2 


2 


Xa 6 \r]\ 2 


(a^aty - A|r?| 2 a 6 )_ 




.(a^aty + \\r]\ 2 a 6 ) 



3 Evolution Of Asymmetry Parameter: 

In order to get some exact results and numerical values for the parameter providing the 
baryon asymmetry of the universe, we consider a (quite realistic) situation where we can 
evaluate the Bogolubov coefficients exactly. 

To this end, we consider the time evolution of wave function under the action of the Hamil- 
tonian Hf given by (42). 
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We go over to the coordinate representation, by defining the operators 



A s (k,t) 
A\(k,t) 
B s (k,t) 
B\(k,t) 



Hf then reduces to the following: 
d 3 k 1 



dt 



n+(t) a° 

n 3 



U A (k,t)+iPu A (k,t)), 



■ r S2 +W^ 



a" 



U A (k,t)-iP UA {k,t)), 



n_(t) 



dt 



(^-^-U B (k,t)+iPu B (k,t)), 



. f o_ (t) 

e — = C-^-u A (k,t)-iP UB (k,t)). 



H f (t) 



(2tt) 3 2 



^) 2 u\(k,t) + p UA (k,t) + (^) 2 n|(M) + Pu B (k,t) 



where 



{ a 6 ' 
(-) 2 



2 _ (^<t> + ^v) 2 



4a 6 



4a 6 



(4a(t) - 1), 
(4/3(0 " !)■ 



(52) 
(53) 
(54) 
(55) 



(56) 

(57) 
(58) 



The time evolution of a wave function x(t) under the action of a Hamiltonian H(t) is 
simply 

H(t) X (t) = ij t x(t), (59) 

From the form of Hf given above, it is clear that it is the direct sum of two independent 
Hamiltonian H A (t) and H B (t) for each of the A s and B s modes. Therefore, the wave function 
x(t) for the Hamiltonian Hf is just the sum of two wavefunctions XA(f) and XB(f) where 
they evolve independently. 

Fourier decomposing the wave function XA(t) and %#(£) and the Hamiltonian H A (t) and 
H B {t), we write 



11 



XAit) = 
XB(t) = 

and each mode evolves as 



d 3 k 



X A(k,t); H A (t) 



(2tt) 3 ' 
d 3 k 

j—^ X B(k,t); H B (t) 



d 3 k 



H A (k,t) 



(2tt) 3 ' 



H A {k,t)xA{k,t) = i—x A (k,t), 
H B (k,t)xB(k,t) =i^ X B(k,t). 



(60) 
(61) 

(62) 
(63) 



Since the Hamiltonians H A and Hb are time dependent harmonic oscillators, in the coor- 
dinate space representation (U A ,P A ) and (Hb, Pb) the wave functions XA(k,t) and XB(k,t) 
can be represented by gaussian.wave functions. 

XA (k,t) = L A (t)e^ w ^ X A(k,0), 
XB(k,t)=L B (t)e^ w ^XB(k,0), 



A time derivative gives 



i^- f XA(k,t) = [ i-JL-itf A W A (k,t) J X A(k,t), 



similarly for XB{k,t) 

while the evolution equation (62) gives 

■ 9 „ , 1 
^ (M) = 2^ 

A similar equation holds for XB(k, t). 
This gives us 



d 2 



XA(k,t), 



W A (t) = -ia 6 /2 



3/nXA(k,t) 



XA(k,tY 



implying 



a n^ A 



X a(K t) + 3-Xa(K t) + ^ X A(k, t) = 0. 



(64) 
(65) 



(66) 



(67) 



(68) 
(69) 

(70) 
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The equations satisfied by the two wave functions for each mode are 

X A (k,t) + 3-x A (k,t) + ^± X A(k,t) = 0, (71) 
a a° 

XB (k,t) + 3-x B (k,t) + ^f X B(k,t) = 0. (72) 

Our aim in this paper has been to calculate the baryon asymmetry of the universe for 
the specific model we have chosen. To this end, we have now obtained the time evolution 
equations for the modes of the wave functions of our baryons and anti-baryons. To solve 
these equations, we shall choose the universe to be FRW as mentioned earlier. 
The metric of our FRW is chosen to be conformally flat. 

ds 2 = a(r) 2 (dr 2 - dx 2 ). (73) 

by employing the scaled time 

dr = a(t) _1 d*. (74) 

The equations of motion for the wave functions X A(k,t) and X B(k,t) given above can be 
transformed into ones that resemble a harmonic oscillator with time dependent frequencies. 

1 ^ ^^4 x#)r)=0) (75) 



a 2 (r) dr a 3 (r) dr dr a 6 

1 d 2 XB(k,r) 2 da dxB{%r) , , 

~TF\ 1 + TT~\1 a + —XB{k, t) = 0. 76 

cl z (t) dr a d (r)dT dr a 3 

A change of variable 

ua = axA (77) 
and similarly for XB(k,r) transforms (75), (76) to 

-u^ + (E-V 1 ( V r))u A = 0, (78) 

and 

-u% + (E + V 2 (rrr))u B = 0, (79) 
where prime denotes differential with respect to r and 



2 



Vi(w) = --^jpi - WlW (81) 
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Mvr) = -irS - [A\ V \*\. (82) 



a(r) c/r 2 

where A = A ^^^ , and u^, are given by the equation (12) and (13). 



These forms of the equations suggest interpretations in two ways. Treating r as a spatial 
variable, they are Schroedinger like equations with E=a;^(r). This allows the calculation of 
reflection and transmission coefficients over potential barriers provided by V\{t]t) and V^t^t) 
terms. On the other hand, they are equations for time dependent harmonic oscillators with 
time dependent frequencies (parametric oscillators). 

In particle creation problems, the potential barrier reflection (R) and transmission (T) 
coefficients can be related to the squeezing parameter (r) through sinh 2 (r) = \v\ 2 = ^. We 
shall use this relationship for the calculation of the squeezing-parameter-dependent number 
operator N(k) = \v\ 2 involved in the evolution of wave functions for the particles and anti- 
particles in our FRW background. 

Following previous work 9\ we assume an oscillating inflaton background 

r\ = Asm(mr), (83) 

where A = |A|e 4 ^ is complex and m is the inflaton mass. Substituting the value of r\ in eq(78) 
and (79) we get the following differential equations, 



u" A + (-— + E - f + ^cos(2mr))u A = 0, (84) 



and 



o!' B B 

u" A + ( + E + — - -cos(2mT))u A = 0, (85) 

where B = \A\ 2 A, 

These are exact equations for the evolution of the baryon and anti-baryon fields in an ex- 
panding universe. Since the number of baryons and anti-baryons are related to the reflection 
and transmission coefficients we see from the above that there will be an asymmetry, since 
the baryons encounter a potential barrier whereas the anti-baryons encounter a potential 
well. 

To analyze these equations for the baryon asymmetry, we first consider the case of constant 
expansion = 0). 

We see then that equations (84) and (85) are Mathieu equations 

B B 

u" A +{E-- + -cos{2mr))u A = 0, (86) 
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and 

B B 

u" A + (E + - - -cos(2mT))u A = 0, (87) 
From the theory of Mathieu equations and parametric resonance we see that the frequency 



2k 2 + rrA + mi + 2 J k 2 + m\J k 2 + m\ 
<4 = * * = (|c) 2 , (88) 

must be half integer multiples of a lowest frequency cu. The resonance for the lowest frequency 
occurs for 

u 2 k - (^) 2 = A n , (89) 

where A n is the width of the the instability band of Mathieu equations. For the values of oo^ 
in this instability band, we assume a broad band resonance such that the Mathieu equation 
has instability bands within which parametric resonance can occur. We shall select the first 
instability region as a broad resonance band. The frequency in a broad resonance band is 
slowly varying, thus 



Since u oc - this implies that 



I — I 



a' 



<cu. (90) 



- < u. (91) 



a 



Thus — being put to zero is a justifiable approximation. 
Now we make another approximation. 

We will assume that m^,m^ u 2 . Thus, to a good approximation, 

uj oc k (92) 

and we find 

u" A + (E - Bsin 2 (mT))u A = 0, (93) 

and 

u" B + {E + Bsin 2 {mr))u B = 0, (94) 

where 

E = k 2 , A = X, B = A|A| 2 (95) 
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In the region of broad resonance we replace the oscillating potential near its zeros with an 
asymptotically flat potential of the form 



\r]\ 2 = \A\ 2 sin 2 (mT) ~ 2\A\ 2 tanh 2 (m [T T '' ). (96) 
In each instability region , we have 

u" A + (k 2 - 2\\A\Hanh 2 { m{ - T ~ Ti) ))u A = 0, (97) 

and 

u'b + + 2X\A\ 2 tanh 2 { m(yT ~ Ti ^ ))u B = 0. (98) 

V2 

We will calculate the transmission and reflection coefficients of the above equations, which 
will then provide the amount of particle production for each case. 

The differential equations (97) can be solved through the following definitions 



k 2 - (2A|A| 2 ) 2 _ m 2 
~p~ 2 ' p ~T 



«i = ^P^Z, p 2 = ^- (99) 

along with a change of variable given by: 

y = p(r- n ) } (100) 

we get the following differential equation 

Ml = 0. (101) 



■ U + <**!£Lw (v ) 



+ 



The transmission |2Z] for this barrier is calculable and in terms of physically known values 



is 



T \=- , -jr • ( 102 ) 



The amount of particle production is given as 
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Now for the anti-particle modes the differential equations (98) governing the evolution of 
their wave functions can also be solved exactly. Defining 




(104) 



and a change of variable as before, 



y = p(j- Ti ) 



(105) 



gives the following reduced differential equation 



cPu 2 
dy 2 



+ 4 + 



(2A|A| 2 ) 
P 2 



tanh 2 (y) u 2 = 0. 



(106) 



The amount of particle production for B s (k, t) modes can be calculated from the transmission 
coefficient and is given by 



We plot the co-wave number dependent number of baryons and anti-baryons in figure 1. 
and figure 2. for various values of the parameter . As expected, from the differential 

equations satisfied by the mode functions, we see that because the baryons see a potential 
barrier and the anti-baryons see a potential well, there is a differential amplification of particle 
and anti-particle modes leading to an asymmetry which survives with time. Therefore, 
from our model we have generated a baryon asymmetry in an entirely non-perturbative 
fashion. The approximations made are only at the end of derivations and have been used to 
illustrate the methodology. It is to be noticed that the equations (84-85) are exact for an 
oscillating inflaton field in an FRW Universe. It is entirely possible to solve these equations 
to incorporate the effects of the expansion of the Universe on the generation of baryon 
asymmetry. Since we are relying on parametric resonance for enhancement /suppression 
of the particle/anti-particle creation, the restriction to the lowest instability band of the 
Mathieu equation is a reasonable one. The total integrated baryon asymmetry of the universe 
is given by 




(107) 




(108) 
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0.1 0.2 0.3 0.4 0.5 0.6 



kbar 



Figure 1: shows the variation of particles (dashed line) and anti particles (solid lines) for 



2A|A| 2 



= .5 as a function of comoving wave number kbar = -. 



. 5 



0.4 



0.3 



0.2 



. 1 



0.1 0.2 0.3 0.4 



kbar 



Figure 2: shows the variation of particles (dashed line) and anti particles (solid lines) for 



2AIAI 



■r- = .25 as a function of comoving wave number kbar = -. 

1 & a 
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4 Conclusions. 

Sakharov's requirements for the generation of a baryon asymmetry in the Universe have 
been the starting point of our paper. Adapting the model of ^3] , to the currently prevalent 
models of inflation and preheating through a background clasical complex scalar inflaton field 
[T7j . we have found equations for the evolution of baryon/anti-baryon modes in an expanding 
FRW metric. We have been able to derive the general evolution equations in an entirely non- 
perturbative fashion. This has been possible through an analogy with quantum optics and 
its methods. In particular, we have found that at any particular time, the Hamiltonian for 
the quantum modes has su(2) and su(l, 1) symmetries. Quantum optical techniques allowed 
us to make these symmetries explicit and thereafter to diagonalize the Hamiltonian. The 
diagonalization procedure further allowed us to relate the rotated and squeezed operators to 
the traditional Bogolubov transformations between the mode function at an earlier time to 
those at later times. The approximations that we have utilized at the end to find explicit 
variation of the number of particles/anti-particles created as a function of the co-moving wave 
number are only therefore illustrative. It is entirely possible to carry out the calculations 
without these approximations. Of course, we will then not have the luxury of the vast 
literature on the solutions of the Mathieu equations and instability regions allowing for 
parametric resonant solutions will have to be looked for. On the plus side, the full equations 
(84-85) allow the possibility of including the contributions to the baryon asymmetry coming 
from the expansion of the Universe. The centerpiece of our effort are the equations [78] and 
[79] , which are exact and show the role played by the expansion factor and the inflaton field 
on the evolution of the particle and antiparticle modes. These can be solved for any type of 
inflationary scenario and will be the subject of future investigations. 
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